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Abstract 

The kaon B parameter is calculated in quenched lattice QCD with the 
Wilson quark action. The mixing problem of the As = 2 four-quark opera- 
tors is solved non-perturbatively with full use of chiral Ward identities, and 
this method enables us to construct the weak four-quark operators exhibiting 
good chiral behavior. We find ^(NDR, 2GeV) = 0.562(64) in the contin- 
uum limit, which agrees with the value obtained with the Kogut-Susskind 
quark action. 



Reliable knowledge of the K — K° transition matrix element Bk is an in- 
dispensable ingredient for further advancement in CP violation phenomenology, 
and the evaluation of this Bk parameter has been one of the major targets for 
lattice QCD calculations. Indeed, lots of effort has been expended towards this 
end. The successful calculations of B K so far achieved exclusively employ 
the Kogut-Susskind quark action that respects chiral U(l) symmetry. Whereas 
the verification whether both Wilson and Kogut-Susskind quark actions yield the 
identical result for Bk is an important step to give a full credit to the lattice QCD 
calculation, the attempts made with the Wilson quark action have not yielded 
much success^, |], |5|: the use of the Wilson action with explicit breaking of chiral 
symmetry causes mixing among four-quark operators of different chiral structure, 
and ensuring the correct chiral behavior of the As = 2 operators is substantially 
more complicated than with the Kogut-Susskind action. Early studies have shown 
that the mixing problem is not adequately treated by perturbation theory, leading 
to an "incorrect answer" for the matrix element || . Attempts were then made to 
solve the mixing problem non-perturbatively with the aid of chiral perturbation 
theory [|]]. Unfortunately, they were not successful since the calculation contains 
large systematic uncertainties arising from higher order effects that survive even 
in the continuum limit. More recently a proposal is made [|[] to improve the chiral 
behavior of the As = 2 operator with the use of non-perturbative renormaliza- 
tion (NPR) ||. The mixing coefficients derived in this approach, however, depend 
sharply on the momentum scale where the operator is to be evaluated, and this 
causes a subtlety in estimating the matrix element. 

In this Letter we propose a non-perturbative method to solve the operator 
mixing problem with use of chiral Ward identities 0. This method fully incor- 
porates the chiral properties of the Wilson action, yielding the As = 2 operator 
that shows good chiral behavior. No effective theories are invoked to estimate 
the matrix element. The final result we obtained is encouraging: it shows good 
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agreement with the result with the Kogut-Susskind quark action. We shall also 
revisit the perturbative method. 

Let us consider a set of weak operators in the continuum {0{\ which closes 
under chiral rotation 5 a Oi = ic^Oj. These operators are given by linear combina- 
tions of a set of lattice operators {O a }, as Oj = J2 a Zi a O a . We choose the mixing 
coefficients Zi a such that the Green functions of {Oi} with quarks in the external 
states satisfy the chiral Ward identity to 0(a). This identity can be derived in a 
standard manner and takes the form 

-2 P z A (J2P a (x)d l (o)I[^(p k )) + 4(0,(0) n^(p.)) (i) 

x k k 

-* £<O<(0) II 4>(Pk)S a 4>(pi)) + 0(a) = 0, 

l k^l 

where pk is the momentum of the external quark, Za and p = (m — 8m)/ are 
constants to be determined from the Ward identities for the axial vector current s[§], 
and P a is the pseudoscalar density of flavor a. We note that the first term comes 
from the chiral variation of the Wilson quark action and the third represents the 
chiral rotation of the external fields. 

The four-quark operator relevant for B K is given by Ovv+aa = VV+AA where 
V = sj^d and A = s^^d. Then, Ovv+aa = VV + AA and Ova = V A form a 
minimal set of the operators that closes under A 3 = diag(l, —1, 0) chiral rotation. 
Taking account of CPS symmetry(note that we take rrid = m s in this article) 0, 
mixing of these operators is written Ovv+aa/^=Zvv+aa (Oq + z\0\ + • • ■ + £404) 
and Ova = Z V a • Z5O5, where the six lattice operators 0, are given by O = 
(VV + AA) /2, O x = (SS + TT + PP) /2, 2 = (SS - TT/3 + PP) /2, 3 = 
(VV - AA) /2 + (SS - PP), 4 = (VV - AA) /2 - (SS - PP) and 5 = VA 
with S = sd, P = sj 5 d and T = s[7 M , j v ]d/2 (/i < v). These operators are in the 
Fierz eigenbasis which we find convenient when taking fermion contractions for 
evaluating the Green functions in (1). 

We consider the four external quarks having an equal momentum p, and denote 
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by Tyy+AA and Tva the sum of the Green functions on the left hand side of 
(HD with external quark legs amputated. Using the projection operator Pj for 
the Fierz eigenbasis corresponding to Oi, we write Tvv+aa/ Zvv+aa = Ts-Ps and 

Tva/Zva = To-Po + ri-PiH hr 4 P 4 . Writing Ovv+aa,va in (0) in terms of lattice 

operators, we obtain six equations for the five coefficients Z\, • ■ ■ , z$: 

T i = ci + c\z 1 + --- + ciz 5 = 0(a), z = (),•••, 5 (2) 

This gives an overconstrained set of equations, and we may choose any five equa- 
tions to exactly vanish to solve for zf. the remaining equation should automatically 
be satisfied to 0(a). We choose four equations to be those for % — 1, - • • , 4, since 
Oi, ■ ■ ■ , O4 do not appear in the continuum. The choice of the fifth equation, i — 
or 5, is more arbitrary. We have confirmed that either To = or Ts = leads to a 
consistent result to 0(a) for z\, ■ ■ ■ , 24 in the region pa<>l. In the present analysis 
we choose T 5 = 0. The overall factor Z V v+aa is determined by the NPR method0. 
We convert the matrix elements on the lattice into those of the MS scheme in the 
continuum using naive dimensional regularization (NDR) renormalized at the scale 
y. = 2GeV§: 

Bx( NDR, „) = ( 1 + ?M ( _41„ g (H) _ 11 + 8 l 0g 2 ) ) g%Mg> (3) 

where p denotes the momentum at which the mixing coeffecients are evaluated. 

For comparative purpose we also calculate Bk with perturbative mixing coef- 
ficients, for which we use the one-loop expression in Ref. || after applying a finite 
correction in conversion to the NDR scheme together with the tadpole improve- 
ment with ajjg(l/a). 

Let us remark here that the equations obtained in the NPR method corre- 
sponds to Tj = for % = 1, • • • ,4 in which the contributions of the first and the 
third term in the Ward identity (1) are dropped, i.e., the full chiral properties are 
not taken into account in the NPR approach. 
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Our calculations are made with the Wilson quark action and the plaquette 
action at /3 — 5.9 — 6.5 in quenched QCD. Table [3] summarizes our run parameters. 
Gauge configurations are generated with the 5- hit pseudo heat-bath algorithm with 
2000(/? = 5.9 and 6.1), 5000(/? = 6.3) or 8000(/? = 6.5) sweep intervals apart. The 
physical size of lattice is chosen to be approximately constant at La « 2.4fm where 
the lattice spacing is determined from m p = 770MeV. At each j3 four values of the 
hopping parameter are adopted such that the physical point for the K meson can 
be interpolated. We estimate m s a/2 from m^/m p = 0.648 for degenerate d and 
s quark masses. Errors are estimated by the single elimination jackknife method 
for all measured quantities throughout this work. 

Our calculations are carried out in two steps. We first calculate Zi and Z vv +aa 
using the quark Green functions having finite space-time momenta. For this pur- 
pose quark propagators are solved in the Landau gauge for the point source lo- 
cated at the origin with the periodic boundary condition imposed on the lattice. 
We next extract B K from the ratio (K°{t = T)d vv+AA (t r )K°(t = l))/{K°(t = 
T)A(t'))/{A(t')K°(t = 1)), each Green function projected onto the zero spatial 
momentum, by fitting a plateau seen as a function of t' . For this calculation quark 
propagators are solved without gauge fixing employing the wall source placed at 
the edge where the Dirichlet boundary condition is imposed in the time direction. 
We obtain at m s /2 by quadratically interpolating the data at the four hopping 
parameters. 

We plot in Fig. 1 a typical result for the mixing coefficients as a function of 
the external quark momenta. The plot shows, as desired, only weak dependence 
of Zi on momentum in the range 0.1^p 2 a 2 ^1.0. This enables us to evaluate the 
mixing coefficients with small errors at the scale p* ~ 2GeV, which always falls 
within the range of a plateau for our runs at (5 = 5.9 — 6.5. We remark that this 
weak scale dependence contrasts to the strong scale dependence obtained with the 
NPR method [|]. 
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In Fig. 2 we compare the mixing coefficients evaluated at the scale p* (filled 
symbols) with the perturbative values obtained with ajjg(l/a) (open symbols) as 
a function of lattice spacing. We remark that a large value of Z2 determined by 
the Ward identities sharply contrasts with the one- loop perturbative result Z2 = 0. 
For the other coefficients, the perturbative results agree with the non-perturbative 
ones in sign and rough orders of magnitude. They differ in quantitative details, 
however. 

Let us examine the chiral property of the operator Oyv+AA by calculating the 
ratio {K°\Ovv+aa\K°) /\(0\P\K°)\ 2 , which vanishes at m q = in the continuum. 
In Fig. 3 we show the results at m q = obtained by a quadratic extrapolation of 
data in m q = (1/K — l/K c )/2, where WI stands for our method using chiral Ward 
identities and PT for tadpole-improved one-loop perturbation theory (numbers are 
given in Table 0). The pseudoscalar density P in the denominator is renormal- 
ized perturbatively for both cases. A significant improvement is clearly seen with 
use of the Ward identities, the ratio becoming consistent with zero at the lattice 
spacing m p a^0.3(a^0.08fm). In the perturbative approach, the chiral behavior 
is recovered only after extrapolation to the continuum limit, where we adopted a 
linear dependence on a expected for the Wilson quark action in the extrapolation 
shown in Fig. 3. 

Our final results for B^-(NDR, 2GeV) are presented in Fig. 4 as a function 
of lattice spacing (see Table 2 for numerical details). The method based on the 
Ward identity (WI) gives a value convergent from a lattice spacing of m p a « 
0.3 (Bk ~ 0.6 — 0.8). The large error, however, hinders us from making an 
extrapolation to the continuum limit. Since the origin of the large error is traced 
to that of the mixing coefficients, we develop an alternative method, which we 
refer to as WIvs, i n which the denominator of the ratio for extracting Bk is 
estimated with the vacuum saturation of Ovv+aa constructed by the WI method. 
In this case the fluctuations in the numerator are largely canceled by those in the 
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denominator, and the resulting error in Bk is substantially reduced as apparent 
in Fig. 4. The cost is that the correct chiral behavior of the denominator is not 
respected at a finite lattice spacing due to the contributions of the pseudoscalar 
matrix element besides the axial vector one. While WI and WIvs methods give 
different results at a finite lattice spacing, the discrepancy is expected to vanish 
in the continuum limit. A linear extrapolation in a of the WIvs results yields 
i?i^(NDR, 2GeV) = 0.562(64), which we take as the best value in the present 
work. This value is consistent with a recent JLQCD result with the Kogut-Susskind 
action, ^(NDR, 2GeV) = 0.587(7) (17) §. 

Intriguing in Fig. 4 is that the perturbative calculation (PT), which gives com- 
pletely "wrong value" at a ^ 0, also yields the correct result for Bk, when extrap- 
olated to the continuum limit a = 0. This is a long extrapolation from negative to 
positive, but the linearly extrapolated value -Br-(NDR, 2GeV)=0. 639(76) is con- 
sistent with those obtained with the WI or Wlys method. We note that a long 
extrapolation may bring a large error in the extrapolated value. 

Finally we mention possible sources of systematic errors in our results from 
quenching effects and uncertainties for Gribov copies in the Landau gauge. With 
the Kogut-Susskind quark action it has been observed that the error due to 
quenched approximation is small |10[ . Whether this is supported by calculations 
with Wilson action we must defer to future studies. For the Gribov problem we 
only quote an earlier study fly] which suggests that ambiguities in the choice of 
the Gribov copies induce only small uncertainties comparable to typical statistical 
errors in current numerical simulations. 

In conclusion our analysis for Bk demonstrate the effectiveness of the method 
using the chiral Ward identities for constructing the As = 2 operator with the 
correct chiral property. We have shown that both Wilson and Kogut-Susskind 
actions give virtually the identical answer for Bk in their continuum limit. We may 
hope that further improvement of our simulations leads to a precise determination 
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of Bk with the Wilson quark action. The application of this method to Bb is also 
straightforward. 

This work is supported by the Supercomputer Project (No.l) of High Energy- 
Accelerator Research Organization(KEK), and also in part by the Grants-in-Aid 
of the Ministry of Education (Nos. 08640349, 08640350, 08640404, 08740189, 
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Tables 



Table 1: Parameters of simulations. 



p 


o.y 


0.1 


o.o 


0.0 


L A x 1 


24 d x 64 


32 J x 64 


40 J x 96 


48 J x 96 


#conf. 


300 


100 


50 


24 


K 


0.15862 


0.15428 


0.15131 


0.14925 




0.15785 


0.15381 


0.15098 


0.14901 




0.15708 


0.15333 


0.15066 


0.14877 




0.15632 


0.15287 


0.15034 


0.14853 


K c 


0.15986(3) 


0.15502(2) 


0.15182(2) 


0.14946(3) 


a' 1 (GeV) 


1.95(5) 


2.65(11) 


3.41(20) 


4.30(29) 


a Ms{ l l°) 


0.1922 


0.1739 


0.1596 


0.1480 


m s a/2 


0.0294(14) 


0.0198(16) 


0.0144(17) 


0.0107(16) 


p* 2 a 2 


0.9595 


0.5012 


0.2988 


0.2056 



Table 2: {K Q \d vv+AA \K°) /\ (0\P\K°) | 2 in the chiral limit and ^(NDR, 2GeV) 
for WI, WI[VS] and PT methods function of (3. 





(K0\O vv+AA \K0) . _ n 


5 X (NDR, 2GeV) 






WI 


PT 


WI 


WI[VS] 


PT 


5.9 


-0.0200(39) 


-0.0415(8) 


+0.38(6) 


+0.168(20) 


-0.468(14) 


6.1 


-0.0068(55) 


-0.0333(10) 


+0.68(11) 


+0.288(29) 


-0.225(22) 


6.3 


-0.0017(74) 


-0.0240(12) 


+0.69(12) 


+0.342(33) 


-0.000(21) 


6.5 


+0.006(10) 


-0.0188(17) 


+0.72(18) 


+0.360(52) 


+0.156(40) 


a = 




-0.0009(31) 




+0.562(64) 


+0.639(76) 
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Figure Captions 

Fig. 1 Mixing coefficients Zi,---,Z4 plotted as a function of external momentum 
squared (pa) 2 for K = 0.15034 at f3 — 6.3. Vertical line corresponds to 
p* « 2 GeV. 

Fig. 2 Comparison of the mixing coefficients zi, •••,24 evaluated at « 2 GeV 
using the Ward identity (WI; solid symbols) and perturbative (PT; open 
symbols) methods. The coefficients are plotted as a function of m p a. 

Fig. 3 Test of the chiral behavior of (K°\6 vv+A a\K ) /\ (0\P\K°) | 2 at m q = for 
the WI and PT methods. The operators are renormalized at 2 GeV in the 
NDR scheme. For both methods we use the same P perturbatively corrected 
with the tadpole improvement. The solid line is a linear extrapolation to the 
continuum limit. 

Fig. 4 K° — K° matrix element 5k(NDR, 2GeV) plotted as a function of the lattice 
spacing for the WI, WI V s and PT methods. The solid lines show linear 
extrapolations to the continuum limit. 
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